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We make an extensive study of evolution of gravitational perturbations of D- dimensional black 
holes in Gauss-Bonnet theory. There is an instability at higher multi-poles i and large Gauss-Bonnet 
coupUng a for D = 5, 6, which is stabilized at higher D. Although small negative gap of the effective 
potential for scalar type of gravitational perturbations, exists for higher D and whatever a, it does 
not lead to any instability. 



PACS numbers: 04.30.Nk,04.50.+h 



INTRODUCTION 



In recent years, higher dimensional black holes have 
been in the focus of high energy physics research. They 
are essential for our understanding of key moments of 
string theory, quantum gravity and brane-world scenar- 
ios. The most important classical property of black holes 
is evidently stability: unstable black holes simply can- 
not exist in our world. In addition, when considering 
higher dimensional black holes in the context of anti-de 
Sitter/Conformal Field Theory (AdS/CFT) correspon- 
dence, that is as a dual backgrounds rather than as real 
black holes, instability means a phase transition in the 
dual theory [l| being essential for understanding the field 
theory at finite temperature and high coupling. An op- 
portunity of creating of mini black holes at particle colli- 
sions in Large Hadron Collider, according to Tev gravity 
extra dimensional scenarios, also gave a strong impetus 
to extensive studies of different properties of higher di- 
mensional black holes, including quasinormal modes 
and Hawking radiation 3. 

The study of evolution of gravitational perturbations 
of the D-dimensional black holes in pure Einstein the- 
ory started from the work of Ishibashi and Kodama [2] 
who managed to reduce the cumbersome perturbation 
equations to the usual Schrodinger-like form. Yet the 
resulting effective potential of the wave equation is not 
always positive definite, so that stability is not taken for 
granted. In a few cases the technique of the so-called 
S-deformations, i.e. of deformations of the wave equa- 
tion which does not touch stability properties, helped to 
transform a potential to a positive definite form, thereby 
proving the stability. In a general case, an extensive nu- 
merical investigation of evolution of gravitational pertur- 
bations allowed to prove stability for higher dimensional 
Reissner-Nordstrom-de Sitter black holes with arbitrary 
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charge and A-term [5[ . 

At the same time, the quantum gravity corrections to 
classical general relativity implies the existence of the 
so-called Gauss-Bonnet term in the dominant order cor- 
rection, that is, the term in the Lagrangian squared in 
curvature. This term vanishes when = 4. Therefore, 
the black holes in the Einstein-Gauss-Bonnet theory Q 
have attracted considerable interest recent years 0. In 
particular, scalar field quasinormal modes of asymptoti- 
cally flat Gauss-Bonnet black holes [9] and of asymptoti- 
cally dS/AdS Gauss-Bonnet black holes were considered 
in [10]. The scalar quasinormal modes of Gauss-Bonnet 
black holes in the regime of high damping were consid- 
ered in The scalar field, propagating in the back- 
ground of a black hole, although gives the qualitative pic- 
ture of evolution of perturbations, is not responsible for 
stability, so that perturbations of Einstein-Gauss-Bonnet 
equations must be considered instead [l^, [13]. Thus, 
Dotti and Gleiser reduced the Einstein-Gauss-Bonnet 
perturbed equations to a wave like form with some effec- 
tive potentials. They found an instability for the two par- 
ticular cases: scalar type of gravitational perturbations 
for I) = 5, and tensor type of gravitational perturbations 
for D = 6, leaving the analysis of stability of general D 
open. The problem was that analytical treatment of sta- 
bility is difficult in many cases: even if we have a wave 
like equation with a potential (extremely cumbersome in 
the Gauss-Bonnet case), negative gaps in potentials can- 
not be easily removed by the S-deformations, because one 
needs to know an ansatz that transforms the potential to 
a positive definite one. 

The aim of our work is to perform a complete numer- 
ical analysis of the evolution of gravitational perturba- 
tions for D-dimensional Gauss-Bonnet black holes with 
D = 5 — 11, what is motivated by string theory and 
quantum gravity, and to determine the stability and in- 
stability regions for these black holes. 

The paper is organized as follows: Sec. II introduces 
the metric and wave like equations for Gauss-Bonnet 
(GB) black holes. Sec. Ill describes the methods of 
time domain integration used here and shows all numer- 
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ical data for evolution of perturbations and quasinormal modes. Sec IV discusses the results obtained. 

I 

II. THE PERTURBATION WAVE EQUATIONS FOR GAUSS-BONNET BLACK HOLES 

The Lagrangian of the Einstein-Gauss-Bonnet action is 

/ = J d'^xV^R +a' J d^Xy^iRabcdR'^'''' - ^RcdR"" + R^). (1) 

Here a' is a positive coupling constant. 
The metric has the form, 

ds^ = f{r)dt^ -J^)- (2) 



^^'^ = ' + a(D-3)(Z)-4) - ' = f + (D-2)Vi 

where a = 167^GDa^ 

In order to measure all the quantities in terms of the black hole horizon ro radius we parameterize the black hole 
mass as 

(3) 

As was shown in [l^, the gravitational perturbations of a Gauss-Bonnet black hole can be decoupled from 
their angular part and reduced to the wave-like equation of the form 

S-|^+nr))*(.,r)=0, dr. = ^^. (4) 

with the effective potentials which have very cumbersome form [l3| , • After some algebra, we managed to simplify 
the potentials for the tensor, vector and scalar types of the gravitational perturbation respectively: 



Vv{r) = f{r)^^-J^A{r) + JrD-^A{r)q{r)^ ^ ^ }^ ^ =, (6) 



r2 V A{r)J ^rD-^A{r)q{r) drl 

^drl ^rD-^A{r)q{r) ' 

V (r) = /(^)^(^) (7) 

' 6ir^{D - 3)M(r)2q(r)8(4cq(r) + {D - l)i?(q(r)2 - 1))2 ' ^ ^ 



We used the following dimensionless quantities 

1/1 1 \ /I 



Air) = 



q(r)2 \2 ' D-3J ' \2 D-3 
B{r) = A{rf (l 



D-A) V D-A 

R 



a{D-2,){D-AY 
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U{r) = b{D-lfR^{l^R)-3{D-lfR{{D-l)R^^24c{l^R))q{r)^ 

+2(L) - l)^{24c{D - l)R^ + 168c2(l ^ R) - {D - l)R^{-?> + + 7D(l + R)))q{rf + 
+2(D - lfR{-lMc^ + - 1)(13 + + c(-84 + 44i? + MD{1 + R)))q{rf + 

- l)^(384c^ - 48c(2 + D{'3D - b))R^ + 192c^(-ll + 1) + (-15 + + 
+ - l)i?^(-3(7 + 55i?) + L>(26 + 106i? + 7D{1 + i?))))g(r)'^ + 
+ - lfR{-QAc^{D - 38) + (I^ - 1)(71 + D{1D - 90))i?2 + 
+16c(303 + 255i? + 13L)^(1 + i?) - 2D(73 + 81i?)))(7(r)^ + 
+4(L> - l)2(96c^(-7 + I^) - 8c(L) - 1)(145 - 74D + QD'^)R^ - 
-8c2(9 - llbR + D(-58 - 34i? + IID{1 + i?))) + (I^ - l)i^^(-5(79 + 23i^) + 
+L>(5(57 + AIR) + L>(-81 - 89i? + 1D{1 + i?)))))g(r)^ - 

-A{D - l)2i?(8c^(43 + (72 - UD)D) ^ {D - l)(-63 + D{99 + L)(-49 + bD)))R^ + 
+4c(321 + 465i? + D{121 - 39R + D(-123 - 107R + 17L)(1 + R)))))q{ry + 
+ (L) - l)(128c^(-9 + D){D - 5) + 32c(L> - 1)(246 + D{9 + D{-55 + 8D)))R^ + 
+64c2(L> - 5)(L)2 - 3 + (49 + (L> - 4)L>)i?) - 

-{D- l)i?2(1173 + 565i? + D{-4{997 + 349i?) + D(6(393 + 217i?) + D{-548 - 452R + 45D(1 + R))))))q{r)^ + 
+ (L) - l)i?(-64c2(D - 5)(36 + D{-13 + 3L))) {D - 1)(635 + L>(-1204 + 3D{294 + D(-92 + 9D))))R^ - 
-8c(L) - 5)(63 + 31R + D(127 + 191i? + D{-47 + + (-79 + D)R))))q{rf + 
+2(L> - 5)(64c^(L> - 5){D - 3) + 8c(L> - l)(-27 + D{U1 + (-43 + D)D))R^ + 
+8c^(D - 5)(-3 + 77i? + D{D -2^{D- IS)R)) ^ {D - lfR^{-33{R - 7) + 
+L)(59 + 43R + L)(-59 - 35R + 9D(1 + R)))))q{ry^ - 

-2{D - 5)R{24c^{-ll + D){D - 5){D - 3) + (L> - if {-65 + L)(81 + L>(7L> - 39)))i?^ + 
+12c(-7 + D){D - 5){D - 3){D - 1)(1 + R))q{ry^ + 

+(L> - 5)2(-l + L>)i?2(16c(26 -^{D- 9)D) ^ {D - 1)(77 - 3i? + L>(-18 + L> + (L> - 2)i?)))g(r)^^ + 
+ (I) - bf{D - 3f{D - lfR^q{rY^, 

A = (I) — 2)(c + 1) = ^(^ + 1^-3) is the eigenvalue of the angular part of the Laplacian. 



III. THE EVOLUTION OF PERTURBATIONS IN TIME DOMAIN 



We study the ringing of GB black hole using a numerical characteristic integration method [is'], that uses the 
light-cone variables u = t — and v = t-\-r^. In the characteristic initial value problem, initial data are specified on 
the two null surfaces u = uq and v = vq. The discretization scheme we used, is 

m) = nw) + nE) - ns) - ^^vmmiiimm + ^(a^) , (s) 

8 

where we have used the following definitions for the points: N = (li + A, + A), W = (ii + A, 'u), E = (ix, + A) and 
S = {u^v). 



To see the correct time-domain profile at late time we 
need to calculate precisely the values of the effective po- 
tential which are used in ([8]). In order to do this we 
must integrate numerically the equation for the tortoise 
coordinate and then solve it with respect to the radial 
coordinate with high accuracy (we worked with a preci- 
sion of ~ 2~^^). We used simple Runge-Kutta method 
for the integration. Since it interpolates the function by 
a cubic spline at each step we are able to find analyti- 
cally r(r^) at each step. The final C++ programm that 
finds the time-domain profiles with arbitrary precision is 



available from the last author upon request. 

Let us note that for D > 4 black holes in ordinary 
Einstein gravity the tensor type of gravitational pertur- 
bations is governed by the same wave equation as test 
scalar field. Therefore both types of perturbations pro- 
duce the same quasinormal mode spectrum, that is they 
are isospectral. This coincidence is remarkable, but does 
not take place for Gauss-Bonnet black holes. Indeed, if 
we compare the tensor quasinormal modes in Table I in 
this paper with test scalar field quasinormal modes of 
0, [13], we can see that quasinormal modes as well as 
effective potentials are quite different. 
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I 


a 


8 


1.346 


10 


1.274 


12 


1.227 


16 


1.170 


20 


1.136 


32 


1.087 


40 


1.070 


50 


1.058 


64 


1.047 




FIG. 1: Threshold a as a function of the inverse multipole number ^. Tensor type of gravitational perturbations D ■ 
The points ^ = 16, 20, 32, 40, 50, 64 were fit by the line a = 2.627-^"^ + 1.005. The theoretical result is at ^ 1.006. 



I 


a 


8 


0.268 


10 


0.258 


12 


0.250 


16 


0.241 


20 


0.235 


32 


0.225 


40 


0.222 


50 


0.219 




FIG. 2: Threshold a as a function of the inverse multipole number I. Scalar type of gravitational perturbations D ■ 
The points I = 16, 20, 32, 40, 50 were fit by the line a = 0.517^"^ + 0.209. The theoretical result is at ^ 0.207. 



TABLE I: Fundamental QNMs for GB black hole perturbation of tensor type {i — 2). a and frequencies are measured in units 
of the horizon radius. 



a 




D = 


= 5 


D = 


= 6 


D = 


= 7 


D = 


= 8 


D = 


= 9 


D 


= 10 


D 


= 11 


0.1 


1 


5898 - 


0.3406i 


2.0112 - 


0.452H 


2.3697 - 


0.554H 


2.7006 - 


0.6492i 


3.0212 - 


0.7342i 


3.3386 


- 0.8092i 


3.6565 


- 0.8753i 


0.2 


1 


6435 - 


0.3289i 


1.9557 - 


0.4269i 


2.2166 - 


0.5246i 


2.4761 - 


0.6105i 


2.7428 - 


0.6827i 


3.0167 


- 0.7429i 


3.2970 


- 0.7936i 


0.3 


1 


6766 - 


0.3157i 


1.8881 - 


0.4086i 


2.0901 - 


0.5034i 


2.3143 - 


0.5820i 


2.5541 - 


0.6456i 


2.8045 


- 0.6958i 


3.0633 


- 0.7340i 


0.4 


1 


6946 - 


0.3014i 


1.8217- 


0.3942i 


1.9874 - 


0.4862i 


2.1907 - 


0.5599i 


2.4126 - 


0.6175i 


2.6448 


- 0.658H 


2.8886 


- 0.6816i 


0.5 


1 


7017- 


0.2870i 


1.7603 - 


0.3823i 


1.9024 - 


0.4718i 


2.0919 - 


0.5423i 


2.2990 - 


0.5953i 


2.5154 


- 0.6242i 


2.7514 


- 0.6299i 


0.6 


1 


7012 - 


0.273H 


1.7046 - 


0.3720i 


1.8308 - 


0.4594i 


2.0100 - 


0.5280i 


2.2031 - 


0.577H 


2.4066 


- 0.5902i 


2.6438 


- 0.5788i 


0.7 


1 


6956 - 


0.2600i 


1.6541 - 


0.3630i 


1.7693 - 


0.4486i 


1.9403 - 


0.5163i 


2.1190 - 


0.5612i 


2.3149 


- 0.5544i 


2.5596 


- 0.5312i 


0.8 


1 


6863 - 


0.2478i 


1.6083 - 


0.3549i 


1.7158 - 


0.4390i 


1.8799 - 


0.5068i 


2.0427 - 


0.546H 


2.2389 


- 0.5176i 


2.4927 


- 0.4894i 


0.9 


1 


6745 - 


0.2366i 


1.5667 - 


0.3476i 


1.6685 - 


0.4303i 


1.8265 - 


0.499H 


1.9716 - 


0.5299i 


2.1767 


- 0.482H 


2.4378 


- 0.4532i 


1.0 


1 


6611 - 


0.2262i 


1.5286 - 


0.3409i 


1.6264 - 


0.4225i 


1.7790 - 


0.493H 


1.9054 - 


0.5095i 


2.1255 


- 0.4497i 


2.3915 


- 0.4219i 
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TABLE II: Fundamental QNMs for GB black hole perturbation of vector type = 2). a and frequencies are measured in units 
of the horizon radius. 



a 


D = 5 






D = 8 


^ = 9 


D = 10 


i:> = 11 


0.1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.8 
0.9 
1.0 


1.0887 - 0.3196i 
1.0490 - 0.3077i 
1.0128 - 0.2950i 
0.9797 - 0.2829i 
0.9494 - 0.2719i 
0.9217 - 0.262H 
0.8961 - 0.2535i 
0.8726 - 0.2458i 
0.8508 - 0.2385i 
0.8302 - 0.2318i 


1.4227 - 0.4372i 
1.3321 - 0.4063i 
1.2582 - 0.3831i 
1.1973 - 0.3654i 
1.1459 - 0.3510i 
1.1014 - 0.3390i 
1.0627 - 0.3287i 
1.0286 - 0.3197i 
0.9982 - 0.3116i 
0.9712 - 0.3042i 


1.7460 - 0.5373i 
1.6016 -0.4964i 
1.4947 - 0.4682i 
1.4113 -0.4466i 
1.3434 -0.4293i 
1.2862 - 0.415H 
1.2369 - 0.4030i 
1.1935 - 0.3927i 
1.1547 - 0.3837i 
1.1194 - 0.3758i 


2.0608 - 0.628H 
1.8664 - 0.57682 
1.7315 - 0.5406i 
1.6287 - 0.51232 
1.5454 - 0.48862 
1.4752 - 0.46762 
1.4143 -0.44792 
1.3606 - 0.42862 
1.3131 - 0.40892 
1.2713 - 0.38862 


2.3712 - 0.70922 
2.1330 -0.64472 
1.9735 - 0.59752 
1.8540 -0.55832 
1.7589 -0.52292 
1.6812 -0.48932 
1.6175 -0.45702 
1.5652 - 0.42662 
1.5220 -0.39872 
1.4859 -0.37352 


2.6805 - 0.78152 
2.4036 - 0.70062 
2.2239 - 0.64012 
2.0923 - 0.58792 
1.9919 - 0.54002 
1.9140 - 0.49652 
1.8525 - 0.45802 
1.8030 -0.42442 
1.7620 - 0.39532 
1.7275 - 0.37002 


2.9919 - 0.84502 
2.6800 - 0.74752 
2.4836 - 0.67272 
2.3448 - 0.60812 
2.2426 - 0.55152 
2.1652 - 0.50292 
2.1045 - 0.46192 
2.0554 - 0.42722 
2.0144 - 0.39782 
1.9794 - 0.37252 



TABLE III: Fundamental QNMs for GB black hole perturbation of scalar type [1 — 2). a and frequencies are measured in 
units of the horizon radius. 



a 










^ = 9 


D = 10 


i:> = 11 


0.1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.8 
0.9 
1.0 


0.8248 - 0.24572 
0.7549 - 0.25662 
instability 
instability 
instability 
instability 
instability 
instability 
instability 
instability 


0.9386 - 0.35572 

? 
? 

0.1906 - 0.29532 
0.2699 - 0.20672 
0.3061 - 0.15152 
0.3264 - 0.11532 
0.3387 - 0.09022 
0.3465 - 0.07222 
0.3513 - 0.05892 


1.1592 - 0.52112 
0.5723 - 0.33762 
0.5939 - 0.21632 
0.6010 - 0.15332 
0.6015 - 0.11642 
0.5989 - 0.09252 
0.5948 - 0.07612 
0.5899 - 0.06422 
0.5847 - 0.05542 
0.5793 - 0.04862 


0.6879 - 0.45922 
0.7443 - 0.27822 
0.7631 - 0.19702 
0.7665 - 0.15252 
0.7639 - 0.12442 
0.7588 - 0.10512 
0.7528 - 0.09102 
0.7463 - 0.08042 
0.7398 - 0.07202 
0.7334 - 0.06532 


0.7935 - 0.42962 
0.8703 - 0.29262 
0.8947 - 0.22812 
0.9018 -0.18932 
0.9020 -0.16292 
0.8990 -0.14352 
0.8946 -0.12852 
0.8895 -0.11662 
0.8840 -0.10692 
0.8785 - 0.09882 


0.8833 - 0.44172 
0.9715 - 0.33822 
1.0044 - 0.28622 
1.0190 - 0.25202 
1.0255 - 0.22682 
1.0278 - 0.20692 
1.0277 - 0.19062 
1.0261 - 0.17702 
1.0236 - 0.16532 
1.0205 - 0.15512 


0.9555 - 0.47502 
1.0501 - 0.40152 
1.0905 - 0.36362 
1.1132 - 0.33742 
1.1277 - 0.31682 
1.1376 - 0.29942 
1.1445 - 0.28412 
1.1494 -0.27042 
1.1527 - 0.25792 
1.1548 - 0.24642 



Here we shall consider uo = corc — iooim-, and the uo 
is chosen so that positive uoim corresponds to a damped 
mode. On the tables I-III one can see fundamental quasi- 
normal modes obtained by time-domain integration. We 
see that imaginary part of uj^ which is proportional to the 
damping rate, is always decreasing, when a is increasing, 
for all three types of gravitational perturbations and all 
D. In other words the less D and the stronger the Gauss- 
Bonnet coupling a, the slower decay of the perturbations 
is. In contrast to the imaginary part, the real oscillation 
frequency ujRe does not behave uniform: uoRe decreases 
as a grows for most cases of tensor and vector modes. 
The behavior of scalar mode is different: here we have 
two competing for the domination modes (see Fig. 6, for 
D = 10) at different stages of the quasinormal ringing. 
Thus for I) = 7, for example, at some values of a the two 
modes with close imaginary parts appear, what makes 
the picture of time evolution more complicated. Then, at 
some larger a the dominance takes another mode. This 
superposition of modes, also with competing excitation 
coefficients of the particular modes, makes dependence 
of the fundamental scalar type QNMs on a and D non- 
monotonic. 

To check that our computation scheme is working 
properly and gives no numerical error, we repeated the 
integration with much smaller step and higher accuracy 
of all incoming data. The picture of the evolution does 
not change, what means high stability and accuracy of 
the integration. Another check may be going to par- 
ticular known limits, such as pure Schwarzschild case 



a = 0. Then we must get pure Schwarzschild QNMs 
in that limit. Again the dominance of the two modes 
at different stages of the ringing complicate the picture. 
Indeed, for instance for 

D = 10 and a = from [5] we have co = 2.45 - 0.98z, 
what is not a fundamental mode for the whole stage 
of ringing but rather for the first period. Indeed, we 
can see the approaching of our GB QNMs to the pure 
Schwarzschild ones in the following data: 



a 






0.00001 
0.0001 
0.001 
0.01 
0.1 


1.2347-0.9328 i 

1.2298-0.93 i 
1.1834-0.9116 i 
0.8960-0.7825 i 
0.8831-0.4418 i 


2.4576-0.9873 i 
2.4580-0.9867 i 
2.4615-0.9807 i 
2.4387-0.9336 i 
2.0611-0.8058 i 



Therefore we can conclude that there is a kind of con- 
ceptual gap of what one can consider as a fundamental 
mode of the quasinormal ringing. Probably, the better 
choice would be to consider the last stage of ringing, im- 
mediately before the tail stage as the one where the fun- 
damental quasinormal modes must be defined. 

Another essential and distinctive feature of the Gauss- 
Bonnet quasinormal ringing is that instability occurs at 
higher multipole numbers while lowest £ are stable! 
This is indeed remarkable as, naively, one would expect 
that if the lowest multipole is stable, then higher mul- 
tipoles just raise up the pick of the potential barrier, so 
that higher multipoles should stabilize the potential. Yet, 
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FIG. 3: Potential and profile for GB black hole perturbation of scalar type (D = 6, / = 2, a = 0.3). The negative gap does not 
lead to instability. It causes exponentially damping "tails" to appear just after the initial outburst. Therefore we are unable 
to see QN ringing. 




FIG. 4: Time-domain profiles for the "region of the irregular QN-ringing" of the gravitational perturbation of scalar type 
D = 6, ^ = 2, a = 0.15,0.20,0.25,0.30,0.35,0.40 (plots from left to right). For a = 0.15 (first plot) we see usual decaying 
oscillations. For a = 0.20 two concurrent modes with the same damping rate. As a increases we observe exponentially damping 
tails, that do not oscillate. At higher a we see the oscillation behavior again, but the frequency of oscillation for a = 0.40 (last 
plot) differs significantly from that for a = 0.15. 



for Gauss-Bonnet black holes, higher multipoles also in- 
crease the negative gap near the black hole horizon [13], 
[l^, allowing existence of bound states in the gap. This 
instability at higher multipoles seems to be intrinsic to 
Gauss-Bonnet theories, because similar instability was 
found some time ago in [17] for Gauss-Bonnet cosmolo- 
gies. 

One can see the evolution of instability in time domain 
in Fig. 5. The larger ^, at the earlier times instability 
growth occurs, and the stronger the growth rate is (Fig. 
5). In addition, the higher the smaller threshold value 
a at which instability happens. At a smaller than the 
critical value, black holes are stable. Therefore it is very 
important not to be limited by small i but to see the 
regime high multipoles, in order to determine the thresh- 



old a with good accuracy. Analytical estimations of [13] 
are compared here with our numerical results in Fig.l, 
2. There one can see that the estimation of p^], which 
gives an expression for the threshold value a ^ Ai~-^-\-B, 
valid for not very large while at larger £ instability oc- 
curs at smaller a. This decreases the minimal value of 
a, when Gauss-Bonnet black holes become unstable. Es- 
sential point is that the above described instability exists 
only for D = 5 and 6, while at higher D the black holes 
are stable. 

It should be explained here that in Table III we have 
two values of a which do not correspond to instabilities, 
yet do not have any definite quasinormal modes. The 
point is that for scalar type of gravitational perturba- 
tions, in addition to the negative gap which may deepen 
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FIG. 5: The picture of instability, developing at large multi- 
pole numbers: D = 6, ^ = 8 (red), i = 12 (green), i = 16 
(blue), a = 1.3. Tensor type of gravitational perturbations. 




and then some other mode start dominating. We stress 
that this rather odd picture, in no way can be suspected 
as an instability, because this "transition" behavior hap- 
pen also for some parameters, for test scalar field, when 
stability can be proved analytically. This kind of nega- 
tive gap at D > 5, is similar to that for the Reissner- 
Nordstrom black holes [5] , which is not deep enough and 
do not produce the instability. 

IV. DISCUSSION 

In this paper the numerical integration in time domain 
was done for the gravitational perturbations of black 
holes in D = 5 — 11 Gauss-Bonnet theories. The instabil- 
ity occurs only for D = 5 and D = 6 cases at some large 
values of a. Higher D stabilize the perturbations. The 
instability starts after some period of quasinormal ring- 
ing and at earlier time for larger multipoles i. Apparently 
the instability at large a is a physically expected result: 
Gauss-Bonnet black holes is inspired by a one-loop string 
theory approximation, so that GB theory is valid only as 
soon as a is small enough. Otherwise one needs to take 
into consideration higher order corrections. 

It is interesting to generalize the present work to the 
case of charged Gauss-Bonnet black holes and asymptot- 
ically de-Sitter black holes, because we know that A-term 
and the black hole charge give also negative gap to the 
effective potential, so that if taking into account all these 
factors, a-coupling, and A and or charge, the parameters 
range of instability might be increased. 



FIG. 6: The picture of time domain evolution for scalar type 
of gravitational perturbations D = 10, ^ = 2, a = 0.01. One 
can see that two modes are dominating at different stages, so 
that in fact the signal is dominated by a superposition of the 
two modes. 

when i is increasing, there is a negative gap at higher 
which is less deep at higher £. This gap does not 
produce instability, but may suppress the quasinormal 
ringing by a exponential tail behavior, so that we almost 
do not see any ringing period (See Fig. 3, 4). For such a 
situation, it would be too strong to state that the appear- 
ing very short period of oscillations with small real fre- 
quency is dominated by some quasinormal modes. This 
"transition" period happens at intermediate values of a. 
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